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Abstract — The next generation of ultra-low-noise cryogenic 
detectors for space science applications require continued ex- 
ploration of materials characteristics at low temperatures. The 
low noise and good energy sensitivity of current Transition Edge 
Sensors (TESs) permits measurements of thermal parameters of 
mesoscopic systems with unprecedented precision. We describe a 
radiometric technique for differential measurements of materials 
characteristics at low temperatures (below about 3 K). The 
technique relies on the very broadband thermal radiation that 
couples between impedance-matched resistors that terminate a 
Nb superconducting microstrip and the power exchanged is 
measured using a TES. The capability of the TES to deliver 
fast, time-resolved thermometry further expands the parameter 
space: for example to investigate time-dependent heat capacity. 
Thermal properties of isolated structures can be measured in 
' geometries that eliminate the need for complicating additional 
' components such as the electrical wires of the thermometer itself. 
Differential measurements allow easy monitoring of temperature 
drifts in the cryogenic environment. The technique is rapid to 
use and easily caUbrated. PreUminary results will be discussed. 



I. Introduction 

The problem of characterizing the thermal properties of 

■ mesoscopic thin-film structures at low temperatures (below 
\ 3 K), in particular their thermal conductances and heat capaci- 

■ ties as a function of temperature, remains one of the key chal- 
I lenges for designers of ultra-low-noise detectors. This problem 
. is not confined exclusively to the detector community. 1 1 1 
"These measurements require small, easily fabricated, easily 

characterized thermometers. Techniques are already used such 
as Johnson noise thermometry (JNT) using thin film resistors 
as noise sources with dc-SQUlD readout or measurements 
of thermal properties using Transition Edge Sensors (TESs) 
but both have practical limitations. JNT can perform measure- 
ments over a reasonable temperature range and is in principal 
a primary thermometer but is in practice secondary because 
of stray resistance in the input circuit to the SQUID that must 
be calibrated. The achievable measurement precision,(TT, for a 



source at temperature Tg is given by the radiometer equation 

T2 



2 

(Trp 



(1) 



where t„i is the measurement time and A/ is the measurement 
bandwidth. 121, Q In practice the bandwidth is limited by the 
source resistance and the input inductance of the SQUID to 
a few lO's of kHz. This gives ax — 3mK for Tg = 500 mK 
with tm — Is. This is the precision that we found in 
practice. f4| TESs can be used to determine conductances by 
measuring the power dissipated in the active region of the 
device where electrothermal feedback (ETF) stabilizes the 
TES at its transition temperature, Tc, as a function of the 
temperature of the heat bath, T^. These measurements are 
widely reported particularly in the context of measurements 
of the thermal conductance of silicon nitride films, but a key 
limitation is that the technique measures thermal properties 
averaged over a large temperature difference (i.e. between Tc 
and Tfc) and the films under study must support additional 
(generally superconducting) metalization to provide electrical 
connection. A problem arises if the thermal properties are 
themselves a function of temperature. A technique for mea- 
suring conductances or heat capacities of micron-scale objects 
rapidly over a reasonable temperature range without additional 
overlying films certainly seems to be required. The proposed 
technique permits true differential measurements of conduc- 
tance (i.e with small temperature gradients), in a geometry 
without complicating additional metalization. Measurements 
of heat capacities are also possible. The technique is easy to 
implement, simple to calibrate and rapid to use. 

We recently demonstrated highly efficient coupling of very 
broadband thermal power between the impedance-matched 
termination resistors of a superconducting microstrip transmis- 
sion line.|5| The efficiency of a short microstrip (I ~ 2mm 
was better than 97% for source temperatures up to 1.5 K. 
When the coupling efficiency is very high the power trans- 
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Fig. L Thermal circuit of a simple geometry. 



ferred along a superconducting microstrip between a source 
and a TES can be used as a fast, accurate thermometer. |6| The 
practical temperature sensitivity is determined by the TES low- 
frequency Noise Equivalent Power, NEP{0). The limiting 
temperature sensitivity is determined by the very-wide radio- 
frequency bandwidth of the microstrip coupling so that we 
substitute A/ Au in Eq.[I]and = 34.5/r^GHz/K is 
the equivalent rf-bandwidth of a blackbody source at temper- 
ature Ts. 

II. TES THERMOMETRY 

We begin by reviewing the technique of TES thermometry 
for a simple geometry. In section III-BI we describe the full 
thermal circuit used in the measurements of conductances. In 
sections III-CI and III-DI we describe how the conductances are 
determined experimentally. 

A. A simple geometry 

A thermal circuit for a simple geometry for TES thermom- 
etry is shown in Fig. [T] A source of total heat capacity Cg 
which could be a thermally isolated silicon nitride (SiN^) 
island is connected to a TES also formed on a SiNx is- 
land by a conductance formed from a resistively-terminated 
superconducting microstrip. Broadband power is transferred 
between the impedance-matched termination resistors of the 
microstrip which has thermal conductance Gm- Over a portion 
of its length the microstrip crosses the Si substrate so that any 
phonon conductance is efficiently heat sunk to the bath. Gm 
arises from conduction due to photons. The source and TES 
are connected to a heat bath at temperature Ti, by conductances 
Gs and Ggb respectively which include the contribution from 
the phonon conductance of dielectric of the microstrip. Heaters 
permit the source temperature to be varied and the TES to be 
calibrated. Changes in the power coupling along the microstrip 
are measured using the TES which is in close proximity to its 
termination resistor For low source temperatures Ts < 3 K, all 
of the power is contained within the pair-breaking threshold 
of the superconducting Nb, 2 A Mb/h 760 GHz where h is 
Planck's constant and A Mb is the superconducting energy gap. 
The power transmitted between the source at temperature Ts 
and the TES at its transition temperature Tc is given by 



where v is the frequency, Pu{Ti) — hvn{v, Ti) and n{v, Ti) is 
the Bose-Einstein distribution, and we have assumed that the 
coupling is loss-less. If the source temperature is low Ts < 3 K 
all of the power is contained well-within the cut-off frequency 
of the microstrip and the upper limit of may be set to infinity. 
Eq. |2] then has the solution 



6h 



(3) 



The measurement conductance Gm{Ts) = dPm{Ts,Tc)/dTs 
is Gm{Ts) — TT^k^Ts/3h. Gm determines how changes in 
source temperature affect power flow to the TES. 

To calibrate the thermometry we need to know the low- 
frequency TES current-to-power responsivity, s/(0). This is 
measured by applying slowly-varying power, SPh, to the 
known heater resistance on the TES island and measuring 
the change in detected current SI. The responsivity is then 
s/(0) = Sl/SPfi- As the source temperature is changed using 
the source heater, the power incident on the TES changes as 
given by Eq. |3] Measuring the change in the current flowing 
through the TES, SI, determines the change in detected power 
SPm = SI/si{0). Noting that 5T,„ = P{Ts,T,) - T(Tb,T,), 
and since ETF fixed the TES's temperature at Tc, we can 
determine the source temperature to a good approximation as 



T, 



' 61 6h 
si{0) TT^kl 



(4) 
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The practical temperature measurement precision is deter- 
mined by the low-frequency TES NEP{0) so that|i6| 

NEP^jO) ^ 

2^mG^(Ts) 

Here we have intentionally omitted thermal fluctuations of the 
source island which we consider a signal in this geometry. For 
NEP{0) = 2 X 10"1^W/\/Hz the achievable temperature 
precision is 30 fiK for i,,, ~ Is and Tg = 500 mK. This is 
two orders of magnitude better than JNT 

B. The measured geometry 

In the full geometry two source islands Si, S2 are connected 
by the subject under test here a conductance, 6*12. The subject 
may be more complicated. Each source island is connected to 
its own TES by a microstrip. Figure |2] shows the full thermal 
circuit. We measure the quiescent TES currents to monitor 
and subtract small drifts in the substrate temperature or the 
electronics. The measurement precision of the temperature of 
either source, ar, with a measurement time tm is determined 
by the low-frequency TES NEP{0) so that 



2 NEP^{0) 



1 



1 



kbTl, 



, (6) 



[P,{Ts) - P,{T,)]di^, (2) 



2tra \GUTb) GUTl,2)J ' Ci,2 

and the temperature measurement precision for this differenc- 
ing approach includes a contribution from the bath temperature 
measurement. We have also now explicitly included the effect 
of thermodynamic fluctuations in the temperatures of the 
sources of heat capacity Ci_2 since these fluctuations directly 
affect the precision with which the source temperature can be 
determined. 
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Fig. 2. Thermal circuit of the full measurement. The subject under test is 
shown here as a conductance G12 but may be more complicated. 
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Fig. 3. Schematic showing the time variation of the input powers (lower 
traces) and output cuiTents (upper traces) for one cycle of the measurement. 
The sequence of measurement and settling times is indicated. The full 
time-sequence represents a measurement frame. The measurement intervals 
determine the temperatures as shown. 



C. Measurement of the conductance G12 

We measure the interconnect conductance G12 in a three 
step process. A schematic of the measurement cycle is shown 
in Fig. [3] In the first step the quiescent TES currents are mea- 
sured. Then dc power Pq is input to both source islands raising 
the temperatures from the bath temperature Tf, to Ti, T2. We 
find Ti = T2 and we quantify the effect of Ti ^ T2 later 
Finally power to ^2 is stepped by an additional small amount 
SP2 raising the island temperatures to T{ and Tj. Since both 
Tl — Ti and — T2 are small the measurement is differential. 
For the next temperature sample the process is repeated with 
incremented Pq. The power steps are adjusted in software 
to give approximately constant increments and differences in 
temperature across the range of temperatures measured. One 
complete measurement cycle (i.e. P = 0, Po, + SP2), with 
a short settling time tgettie between each measurement step to 
accommodate thermal response times, defines a 'frame' time. 
The frame must be measured in a time less than the Allan 
time of the system. 

The power flow across a conductance G connecting thermal 
reservoirs at temperatures T, T' can be written for notational 
convenience as 



P{T',T)^ [ G{T)dT ^ G{T',T){T' -T), 
Jt 



(7) 



where the over-set line denotes averaging. If T' — T = ST 
is small then the power flow can be linearized so that 



P{T', T) = G{T)5T. Ignoring for now the small conductance 
of the microstrips, the input powers and resultant temperatures 
are related by 



-Po - Pib{Ti,Tb) - Gi2{Ti,T2){T2 - Tl) (8a) 

Po = Pib {T2 , Tfc) + Gi2{Ti,T2){T2 - Tl ) (8b) 

Po = Pifc(Ti', Tfc) - Gi2(T{,T^)(T^ - T[) (8c) 

To + m = T2fc(T^, Tfc) + Gi2(T{,T^)(T^ - Ti) (8d) 



Subtracting [8a| from [8cl [8b] from [8d] and using, for example, 

T2b(T^,Tfc)-T2fc(T2,Tfc)) = / G2b{T)dT 

= G2b{%)ST2, 



(9) 



where T2 = +T2)/2 and the final equality follows since 
ST2 is small, we find 

SP2 = 2Gi2iT^){ST2 - STi) + G2bi%)6T2 - Gu{%)5T^ 
_ (10) 
and T12 (Tl + T[ + T2 + T^)/4. Finally the effect of 
conductance along the microstrip needs to be included. The 
result is 

Gi2(Ti2) = 

^T2 - (G26(7^) + G„(7^)) 6T2 + (Gi6(7Y) + G™(7Y)) ST^ 



2(<5T2 - ,5Ti) 



(11) 



D. Measurement of the conductances Gib, G26 

The conductance to the bath at a given temperature is 
measured in a two-step procedure where we measure the 
quiescent TES current (P — 0) and the effect of applying 
equal power Pq to both islands and measuring the changes 
in Tl and T2. The next temperature sample uses incremented 
To- This measurement is rapid. With t„i = 0.82s, 500 data 
points are acquired in less than 15 minutes. In the analysis we 
use Eqs. [8al and [8bl and assume Ti — T2. This introduces an 
inevitable error in the analysis and the magnitude is of order 
e = |Gi2(Ti2)(T2 - Ti)|/Pi,2b. Experimentally the error is 
small ^ 1%. A fifth-order polynomial is fitted to the T1.2 — To 
data and the conductances Gib, G26 found by differentiation. 
Note that since the temperature difference T1.2 — T^ may be 
large this measurement is temperature-averaged. 

III. Measured device 

An optical image of the device described here is shown 
in Fig. m The measured conductance labelled G12 is a long 
thin SiNx bar of dimensions 500 x 10 x 0.5 /xm'^ formed 
by reactive ion and deep reactive ion etching of a nitride- 
coated Si wafer G12 carries no additional films. G12 connects 
two larger 0.5 /xm-thick nitride source islands Si and S2 
themselves isolated from the Si wafer by four supporting 
nitride legs, two of length 255 /im two of length 358 fim each 
of width 15 /im. One of the longer legs also carries a Nb 
microstrip line. The microstrip is terminated by an impedance- 
matched AuCu termination resistor at each end. Each island 
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Fig. 4. Optical image of the device under test. The central nitride bar 
(labelled G12) of length 500 /^m, width 10 /im connects the nitride source 
islands Si, 82- Resistively terminated microstrip lines run from the islands 
to TESs in the upper comers of the image. The conductance between each 
nitride island and the heat bath is formed from four nitiide legs, two of length 
285 fim, two of 355 fim that run diagonally . 



also supports a square AuCu resistor with Nb bias lines that 
can be used as a heater to modulate the temperature of the 
sources. The microstrips run over the Si wafer then onto nitride 
islands which support TESs. Routing of the microstrips in this 
way ensures that the phonon conductance associated with the 
microstrip dielectric is efficiently heat-sunk to the bath. The 
TES islands also include AuCu resistors with Nb bias lines 
that allow the power-to-current responsivity of the TESs to be 
measured. 

AuCu resistors are 40 nm thick, Nb bias lines and the 
ground plane of the microstrip are 250 nm thick. The mi- 
crostrip dielectric is sputtered Si02 and is 400 nm thick. The 
TESs use our standard higher temperature MoCu bilayer layup 
with 40 and 30 nm of Mo and Cu respectively. The fabrication 
route is identical with our usual process for MoCu TESs.|7| 
The TESs are voltage biased and read-out with SQUIDs. 
The device is measured in a He-3 refrigerator with a base 
temperature of 259 mK. The transition temperature of the 
TESs described here was Tc ~ 485 mK which is slightly 
higher than reported in our earUer work with the same Mo-Cu 
bilayer lay-up. As a result the conductance to the bath of the 
TESs is increased and the measured TES Noise Equivalent 
Power is NEP{0) = 1.2 x lO^^^ 

IV. Results and Discussion 

The Allan variance statistic, (7^(t) where t is the inte- 
gration time, provides an exceptionally powerful diagnostic 
of system stability and a plot of the variance as a function 
of r identifies the optimum time for signal averaging. [ 8 1, 
[|9l In the Allan plot, a log-log plot of a\ as a function 
of r, underlying fluctuations with frequency-domain power 
spectra varying as 1//" show a r^""^' dependence. Hence 
white noise with a = has a r^^ characteristic. 1/f noise 
shows no dependence on integration time and drift exhibits 
a dependence with 1 < a < 3. Figure |5] shows measured 
Allan variances for both SQUID systems with biased TESs 
and the bath temperature at Tb = 259 ±0.5 mK. The measured 
characteristic indicates that white noise is reduced by time- 
integration up to a maximum of order r^i ~ 5 s. The dashed 
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Fig. 5. Allan plot for the two SQUID systems with biased TESs and 
Tj, = 259 lb 0.5 mK. The optimum integration time, identified as ta, is at 
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Fig. 6. Difference in temperature between the two islands as a function 
of average temperature for equal applied powers. Note the improvement in 
the measurement precision as the temperature increases. The inset shows the 
measured temperature precision for repeated measurements near 384 mK. The 
measured precision ax = 0.54 mK compares well with the calculated value 
from Eq. |6] which gives or = 0.47 mK for the measurement time used 
im = 0.82s and the measured NEP{0) = 1.2 X IQ-^'^ WHz-l/^ 



lines show the expected behaviour for white noise with a 
slope and for drift with a slope of r^. This shows that 
drift limits these measurements. Guided by the Allan plot, 
we chose a sample time t,n = 0.82 s being 2^^^ data points 
sampled at 20 kHz and tsettie = 100 ms. For the conductance 
measurements with three power steps the total frame time is 
2.8 s at each sample temperature. 

Applying equal power to both sources, Gu and G26 were 
determined. As discussed earlier the analysis introduces an 
unavoidable error in the measurement of Gu, and G2b depend- 
ing on T2 — Ti. Figure |6] shows the difference in temperature 
6T12 ^ T2 — T1 for equal powers applied to both islands as a 
function of average island temperature. The difference is small 
and a maximum of about 5 mK at 500 mK. Since ST12 is small 
we will see that the error is negligible. The temperature differ- 
ence implies a difference in conductance for the two notionally 
identical conductances Gih, G26 of about ±2% between 0.26 
and 1.3 K. The temperature dependence of ST12 evident in 
Fig.|6]is also unexpected. It does not seem that the difference 
can be accounted for by experimental uncertainty (such as 
in the calibration of the TES responsivities). One possibility 
might be differences in the actual coupling efficiencies of the 
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Fig. 7. Measured conductances to the bath Gn,, G2b ™d calculated sample 
conductance G12. Estimated en'ors for G12 are indicated. The dotted hues 
indicate dependencies proportional to and T. The single circle is the 
measured Gtes scaled by the A/l ratios to compare directly to Gift and 
G2b- 



microstrip lines. 

The inset of Fig. |6] shows a histogram of calculated tem- 
peratures for repeated measurements of a source temperature 
near 384 mK. The measured precision <tt — 0.54 mK com- 
pares well with the value calculated from Eq. |6] which gives 
at = 0.47 mK for the measurement time used t„i ~ 0.82 s 
and the measured NEP{0) = 1.2 x lO^^^ WKz^i/^ 

Figure |7] shows measured conductances Gu, G2b and 6*12 
as a function of temperature. Representative error bars for 
G12 determined from the variance of repeated measurements 
are indicated. We can now estimate the magnitude of the 
error in Gu and G2b- The maximum temperature difference 
occurs near 500 mK where 6T12 ~ 5 mK is greatest. We find 
e ~ 0.2%, but less than this over most of the temperature 
range. This is considered acceptable. We also show the vari- 
ation with temperature if G = kT^ with k a constant. The 
dotted lines show dependencies /? = 1 and 2. There is a strong 
suggestion here that a simple power law does not account for 
the conductance across this measurement range. At the highest 
temperatures /3 > 2, at the lowest /3 < 2. A reduction of the 
exponent may be expected at low temperatures if dominant 
phonon wavelengths start to become comparable to the nitride 
thickness. The single dot in Fig. |2]is the measured Gtes for 
one of the TESs obtained in the standard way by measuring 
the power plateau in the ETF region of the current-voltage 
characteristic as a function of the bath temperature. This single 
point represents approximately 2 hours of data taking. By 
contrast the conductances to the bath Gib, G26 plotted in Fig. [7] 
are found from 500 measurements of temperature acquired in 
15 minutes. A comparable acquistion time measures G12 with 
50 data points with additional signal averaging. 

Figure [8] shows the ratio of G12 to Gib as a function of 
temperature. The simplest model for power flow along, or 
the conductance of, a uniform bar would assume that the 
power scales as the ratio of cross-sectional area. A, to length, 
I. From the dimensions of the nitride support legs for the 
source islands and the dimensions of G12 we calculate a 
conductance ratio of 0.082. This ratio includes the conductance 
associated with the Si02 dielectric of the microstrip lines and 
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Fig. 8. Ratio of the measured conductance G12 to the conductance to the 
bath Git- The aiTow indicates the expected ratio if the conductances scale as 
the area to length ratio of the structures. 

the smaller contribution from the Nb wiring all of which 
we assume are equal to that of SiN^. The variation of the 
ratio is experimentally significant and may already illustrate 
the difference between thermal properties averaged over a 
temperature range and the true differential measurement. 

V. Conclusion 

We have described our first true-differential measurements 
of thermal conductances of a micron-scaled object at low 
temperatures using microstrip-coupled TES thermometry. The 
temperature precision is akeady significantly greater than that 
achievable with JNT with the same measurement time. The 
measurements are rapid and easily calibrated. The achieved 
precision akeady strongly suggests that the thermal transport 
characteristics of the nitride structure are not described by a 
simple power-law across the temperature range 0.26 to 1.5 K. 
The device under test can be fabricated without additional 
metalization for wiring. It should be straight-forward to in- 
clude specific layers on the nitride test structure to measure 
particular thermal properties in a controlled manner In the 
future we expect to be able to explore the temperature depen- 
dence of heat capacities of thin films such as Si02 including 
possibly measurements of time dependent heat capacity. We 
will explore the effect on conductance of superposed layers: 
for example, does thermal conductance in thin multilayers 
really scale as total thickness? We also see the possibility of 
exploring the engineering of the nitride to realise phononic 
structures, to achieve reductions of the conductance in compact 
nitride structures needed for the next generation of ultra-low- 
noise detectors. 
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